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We evaluate thermodynamical quantities and the transport coefficients of a dense and 
hot hadronic matter based on the event generator URASiMA (Ultra-Relativistic AA col- 
lision Simulator based on Multiple Scattering Algorithm) with periodic boundary condi- 
tions. As the simplest example of the transport coefficients we investigate the temperature 
dependence and the chemical potential dependence of the baryon diffusion constant of a 
dense and hot hadronic matter. 



1. Introduction 



Physics of a high density and high temperature hadronic matter has been highly stud- 
ied in the context of both high energy nuclear collisions and cosmology JIJ. In recent 
ultra-relativistic nuclear collisions, though the main purpose should be confirmation of 
Quark-Gluon Plasma(QGP) state, physics of hot and/or dense hadronic state dominates 
the system. Hence, thermodynamical properties and transport coefficients of a hadronic 
matter are essentially important understand the space-time evolution of the exited region. 
In the cosmology, physics of hot hadrons is important not only for the global evolution of 
the early universe, but also in the nucleosynthesis problem baryon diffusion would play 
an important roll. 

Because of the highly non-perturbative property of a hot and dense hadronic state, the 
thermodynamical properties and transport coefficients have not been thoroughly investi- 
gated. In this paper, we evaluate the transport coefficients by using statistical ensembles 
generated by Ultra-Relativistic A-A collision simulator based on Multiple Scattering Algo- 
rithm (URASiMA). Originally, the URASiMA is an event generator for the nuclear colli- 
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sion experiments based on the Multi-Chain Model(MCM) of the hadrons||. We improved 
the URASiMA to recover detailed balance at temperature below two hundred MeV, and 
we can obtain natural thermodynamical behavior without Hagedorn-type strange tem- 
perature saturation. This is the first calculation of the transport coefficient of a hot and 
dense hadronic matter based on an event generator 0. 

2. URASiMA for Statistical Ensembles 

The URASiMA is a relativistic event generator based on hadronic multi-chain model, 
which aims at describing nuclear-nuclear collision by the superposition of hadronic col- 
lisions. Hadronic 2-body interactions are fundamental building blocks of interactions in 
the model, and all parameters are so designed to reproduce experimental data of hadron- 
hadron collisions. Originally, the URASiMA contains 2-body process (2 incident particle 
and 2 out-going particles), decay process (1 incident particle and 2 out-going particles), 
resonance (2 incident particles and 1 out-going particle) and production process (2 inci- 
dent particles and n (> 3) out going particles). 

The production process is very important for the description of the multiple production 
at high energies. On the other hand, in the generation of statistical ensembles in equilib- 
rium, detailed balance between processes is essentially important. If re-absorption process 
does not exist, one-way conversion of energy into particle production occurs rather than 
heating up and results in artificial saturation of the temperature 0. Therefore, role of 
re-absorption processes is very important and we should take it into account. However 
exact inclusion of multi-particle re- absorption processes is very difficult. In order to treat 
them effectively, multi-particle productions and absorptions are treated as 2-body pro- 
cesses including resonances with succeeding decays and/or preceding formations of the 
resonances. Here two body decay and formation of resonances are assumed. For example, 
NN — > NNn is described as NN — > NR followed by decay of R — > Nn, where R denotes 
resonance. The reverse process of it is easily taken into account. In this approach, all 
the known inelastic cross-sections for baryon-baryon interactions up to y/s < 3GeV, are 
reproduced. 

For a higher energy, y/s > 3GeV, in order to give appropriate total cross section, we 
need to take direct production process into account. Only at this point, detailed balance 
is broken in our simulation; nevertheless, if temperature is much smaller than 3 GeV, the 
influence is negligibly small. 



Table 1 



Baryons, mesons and their resonances included in the URASiMA. 
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In order to obtain an equilibrium state, we put the system in a box and imposed a periodic 
condition to the URASiMA as the space-like boundary condition. Initial distributions of 
particles were given by uniform random distribution of baryons in a phase space. Total 
energy and baryon number in the box were fixed at initial time and conserved through-out 
the simulation. Though initial particles were only baryons, many mesons were produced 
through interactions. After a thermalization time-period of about 100 fm/c, the system 
seemed to be stationary. In order to confirm the achievement of equilibrium, we calculated 
energy distributions and particle numbers. Slope parameters of energy distribution of all 
particles became the same value in the accuracy of statistics. Thus, we may call this value 
as the temperature of the system. The fact that numbers of species saturate indicates the 
achievement of chemical equilibrium H. 

Running the URASiMA many times with the same total energy and total baryons in 
the box and taking the stationary configuration later than t = 150 fm/c, we obtained 
statistical ensemble with fixed temperature and fixed baryon number (chemical potential). 
By using the ensembles obtained through above mentioned manner, we can evaluate 
thermodynamical quantities and equation of states||. 

3. Diffusion Constant 

According to the Kubo's Linear Response Theory, the correlation of the currents 
stands for admittance of the system(first fluctuation dissipation theorem) and equiva- 
lently, random-force correlation gives impedance(Second fluctuation dissipation theorem) 
As the simplest example, we here focus our discussion on the diffusion constant. 
First fluctuation dissipation theorem tells us that diffusion constant D is given by cur- 
rent (velocity) correlation, 

1 r 00 

D = - < v(t) ■ v(t + t') > dt'. (1) 
3 Jo 

If the correlation decreases exponentially, i.e., < v(t) ■ v(t + t') >oc exp (— -), with r 
being relaxation time, diffusion constant can be rewritten in a simple form, 

D = \ < v(t) ■ v(t) > r. (2) 

Because of the relativistic nature of our system, we should use j3 = — = ^ instead of v 
in eq.(|l|) and D is obtained by, 

1 roo 

- / < /3(f) - /3(f + f') > dt'c 2 , (3) 
3 Jo 

\ < /3(f) • (3(t) > c\ (4) 
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with c being the velocity of light. Figure 1 shows correlation function of the velocity 
of baryons. The figure indicates that exponential damping is very good approximation. 
Figure 2 displays the our results of baryon diffusion constant in a hot and dense hadronic 
matter. 
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Figure 1. Velocity correlation of the baryons as a function of time. Lines correspond to 
the fitted results by exponential function. Normalizations of the data are arbitrary. 
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Figure 2. Diffusion constant of baryons. 
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Our results show clearer dependence on the baryon number density while dependence 
on temperature is mild. This result means that, for the random walk of the baryons 
in our system, baryon-baryon collision process is more important than baryon-meson 
collisions. In the inhomogeneous big-bang nucleosynthesis scenario, baryon-diffusion plays 
an important roll. The leading part of the scenario is played by the difference between 
proton diffusion and neutron diffusion |7[]. In our simulation, strong interaction dominates 
the system and we assume charge independence in the strong interaction, hence, we can 
not discuss difference between proton and neutron. However obtained diffusion constant 
of baryon in our simulation can give some kind of restriction to the diffusion constants of 
both proton and neutron. 



4. Concluding Remarks 

We evaluate diffusion constants of baryons in the hot and dense hadronic matter with 
use of statistical ensembles obtained by an event generator URASiMA. Our results show 
clear dependence on baryon number density and weak dependence on temperature. The 
temperature in our simulation is limited to a small range, i.e., from 100 MeV to 200 MeV, 
and this fact can be one of the reasons for unclear dependence on temperature. Strong 
baryon number density dependence indicates that, for the baryon diffusion process, baryon 
plays more important roll than light mesons. In this sense our simulation corresponds to 
high density region and a non-linear diffusion process occurs. 

Calculation of the diffusion constants is the simplest example of first fluctuation dissi- 
pation theorem and the entrance of the non-equilibrium physics. From diffusion constant, 



1.3 



1.1 



• e = 0.16fm"; 
□ e = 0.23 frrT 

Ae = o.3i fn-r 



0.9 



> 

a o.7 

to 
O 



Oot 



0.5 



0.3 



80 100 



120 140 
T [MeV] 



160 



180 



Figure 3. Baryon number charge conductivity as a function of temperature 
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we can calculate charge conductivity |§. Figure 3 shows baryon number conductivity a&, 

cs = ^D, (6) 

where tib is baryon number density, T is temperature and fee is Boltzmann constant (put 
as unity through out this paper), respectively. Therefore, if we want, we can discuss Joule 
heat and entropy production in the Baryonic circuit with use of the above baryon number 
conductivity. 

In principle, taking correlation of appropriate currents, i.e. energy flow, baryon number 
current, stress-tensor, etc., we can evaluate any kinds of transport coefficients. Though, 
in relativistic transport theory, there exist several delicate points[|], once we establish 
phenomenological equations for the high temperature and high density hadronic matter, 
we can evaluate the appropriate transport coefficients in the same manner. Detailed 
discussion will be reported in our forthcoming paper. 
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